Lattice model of three-dimensional topological singlet superconductor with 

time-reversal symmetry 
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We study topological phases of time-reversal invariant singlet superconductors in three spatial 
dimensions. In these particle-hole symmetric systems the topological phases are characterized by an 
even-numbered winding number v. At a two-dimensional (2D) surface the topological properties of 
this quantum state manifest themselves through the presence of v flavors of gapless Dirac fermion 
surface states, which are robust against localization from random impurities. We construct a tight- 
binding model on the diamond lattice that realizes a topologically nontrivial phase, in which the 
winding number takes the value v — ±2. Disorder corresponds to a (non-localizing) random SU(2) 
gauge potential for the surface Dirac fermions, leading to a power- law density of states p{e) ~ e^''^. 
The bulk effective field theory is proposed to be the (3-1-1) dimensional SU(2) Yang-Mills theory 
with a theta-term at 6* = vr. 

PACS numbers: 73.43.-f, 73.20.At, 74.25. Fy, 73.20.Fz, 03.65.Vf 



Bloch- Wilson band insulators are commonly believed 
to be simple and well understood electronic states of mat- 
ter. However, recent theoretical ||l|, |^, || ^ ^, ||] and ex- 
perimental iQ, ^ progress has shown that band insulators 
can exhibit unusual and conducting boundary modes, 
which are topologically protected, analogous to the edge 
states of the integer quantum Hall effect (QHE). These 
so-called Z2 topological insulators (also known as 'quan- 
tum spin Hair insulators), which exist in two- and three- 
dimensional (3D) time-reversal invariant (TRI) systems, 
are characterized by a topological invariant, similar to 
the Chern number of the integer QHE. Given these newly 
discovered topological states, one might wonder whether 
there exists a general organizing principle for topological 
insulators. Indeed, the integer QHE and the Z2 topo- 
logical insulators are in fact part of a larger scheme, dis- 
cussed in Ref. which provides an exhaustive classi- 
fication of topological insulators and superconductors in 
terms of spatial dimension and the presence or absence 
of the two most generic symmetries of the Hamiltonian, 
time-reversal and particle-hole symmetry Q . 

Using this classification scheme which was originally 
introduced in the context of disordered systems [|oj, it 
was shown in Ref. ^ that, besides the 3D Z2 topological 
insulator, there are precisely four more 3D topological 
quantum states. Among these there is one which is par- 
ticularly interesting from the point of view of possible 
experimental realizations. It is called the topological su- 
perconductor in symmetry class CI in the terminology of 
Ref. 1^ and can be realized in time-reversal invariant sin- 
glet BCS superconductors. While the bulk is fully gapped 
in this topological quantum state, there are gapless ro- 
bust Dirac states at the two-dimensional boundary. The 
CI topological superconductor is unique among 3D topo- 
logical quantum states in that it does not break SU(2) 



FIG. 1: (Color online) The diamond lattice has two inter- 
penetrating Bravais sublattices (FCC) colored in blue (light) 
and red (dark), respectively, (a) The pairing amplitude and 
(b) the second neighbour hopping are schematically shown, 
where bonds colored by red (dark) and blue (light) have pair- 
ing potentials or second neighbour hopping amplitudes which 
differ by a sign. 



spin rotation symmetry, and therefore supports the trans- 
port of spin through gapless surface modes. The different 
topological phases of the CI topological superconductor 
can be characterized by an integer winding number v, 
that takes on even values only, and which can be inter- 
preted as the number of species of gapless surface Dirac 
fermions. 

In this Letter, we construct a lattice BCS Hamiltonian 
which realizes both the topologically trivial and the non- 
trivial phases of the CI topological superconductor. We 
compute the winding number and demonstrate the exis- 
tence of gapless two-dimensional Dirac fermions at the 
boundary. The model we consider is defined on the di- 
amond lattice with (the lattice analogue of) spin singlet 
d-wave pairing, and has the form of a 4 x 4 Bogoliubov- 
de Gennes (BdG) Hamiltonian. In reciprocal space the 
noninteracting Hamiltonian reads H = J2k ^fc^(^)^fc' 
with ^E* 
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and bk., 



represent the electron annihilation operators with spin a 



2 



and momentum k on sublattice A and B of the diamond 
lattice, respectively, and 
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term is 8fe = J2i 



Here, the nearest neighbor hopping term is given by 
the next nearest neighbor hopping 
^/_^ik-{si-sj) _|_ g^j^^ ^Yie pairing 

potential is = X^i^^j A^ e'^^'^^'^^j^ where Si=i...4 de- 
notes the four first neighbor bond vectors. The near- 
est and second nearest neighbor hopping amplitudes are 
parametrized by the vector ti and the symmetric matrix 
t'j^j , respectively. Similarly, the symmetric matrix Ay de- 
notes the singlet BCS pairing order parameter, whereas 
Us is the staggered chemical potential. 

From Eq. (|l|), the energy eigenvalues = 
i-v/l^fcP + + lA/ep are readily obtained, exhibit- 
ing a two-fold degeneracy for each k. The symmetry 
operation that realizes particle-hole symmetry (PHS) for 
a singlet pairing BdG Hamiltonian is given by |^ , 

ryn^{-k)ry = -nik), (2a) 

where Vy is the second Pauli matrix acting on the particle- 
hole space. Hamiltonian (^ automatically satisfies sym- 
metry property (|2^), as can be checked easily. If fur- 
thermore time-reversal symmetry (TRS) is present, Ti.{k) 
obeys 



n*{-k)=nik), 



(2b) 
are all 



which is the case, if the pairing amplitudes A^^ 
purely real. The discrete symmetry constraints (|2a 
and (2b) define the CI symmetry class in the Altland- 
Zirnbauer classification It is important to note, 

that an arbitrary Hamiltonian belonging to symmetry 
class CI can be brought into block off-diagonal form. 
This is achieved by means of a unitary transformation 
which rotates the r„, matrices such that (fa, , 7\, , ) — > 

trans- 



forms into 



• zi ' y)- 



Under this rotation TC{k), Eq. 



/ D{k) 
\D^k) 



nik) 

where the upper right block is given by 
Dik)-- 



Afc - iOk -i*fc 



(3a) 



(3b) 



which satisfies D^i—k) = Dik), since A_fc = A^ and 



In order to define a topological invariant for the CI 
topological insulator, we need to introduce, following 
Ref. the projection operator Qik) 



g(fc) = 



qik) 
q^ik) 



qik) = 



E+ik) ' 



(4) 



where, as a consequence of TRI and PHS, q'^i—k) = 
qik). The block off-diagonal form of nik), and hence 
of Qik), is essential to uncover the topological structure 
of the space of all possible quantum ground states in 
class CI. It allows us to introduce a topological invariant 
that classifies maps from the Brillouin zone (BZ) into the 
space of the projection operators (5(k). This invariant is 
a winding number defined in Ref. Q through 

'^'^ e'^-^Ptr[iq-^d,q) (q-'d.q) iq~'d,q)] , 



247r2 



(5) 

where the integral is over the 3D BZ [0. Due to the 
class CI constraint, q^i—k) = qik), v can take on only 
even integer values. (In the absence of any constraint, 
i.e., for symmetry class AIII the winding number can 
be an arbitrary integer.) The winding number changes 
only when the quantum system undergoes a quantum 
phase transition, which is accompanied by the closing of 
the bulk gap. 

We now turn to a more detailed specification of our 
lattice model. When ti = t for alH = 1, . . . , 4 



(6) 



The pairing potential we consider is "(i3^2_j,2"-like, where 
the pairing amplitude on bonds within the x — y plane 
differs in sign form the out-of-plane pairing amplitudes 
[see Fig. 11(a)] 



Afc = 4A 
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For definitiveness we will take t and A positive through- 
out the paper. The BCS dispersion given by and A^ 
has four point nodes (Dirac points), at which the bands 
are four-fold degenerate 

K\± = 27r ( ±1/3, 1, ) , R\± = 27r ( 1, ±1/3, ) . 

Note that PHS and TRI relate Ka^+ to Ka- (a = 1,2), 
as Ka^^ is identified to —Ka.+ via a translation by re- 
ciprocal lattice vectors. The degeneracy at the point 
nodes can be lifted by a finite staggered sublattice po- 
tential or by second neighbor hopping amplitudes, which 
are parametrized by the symmetric matrix t'^y In order 



to open up a bulk gap, we choose t'^ 



13 



— <23 = t' ^ 0, while ~ = 0, in which case the sec- 
ond neighbor hopping term simplifies to [see Fig. |l|(a)] 

kx 



0fc = 4t' cos — I cos — 
2 V 2 



cos — ] + 



(8) 



where we have also included a staggered chemical poten- 
tial. 

The BdG Hamiltonian (|^) has four energy bands, two 
of which arc occupied. With the choice for the hopping 
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FIG. 2; (Color online) (a) Phase diagram as a function of 
second neighbor hopping t' and staggered chemical potential 
/Xs. (b) Numerical evaluation of the winding number. 



parameters and gap amplitudes given by Eqs. the 
spectrum has a bulk energy gap, except when ±6t' = /i^, 
i.e., when the system undergoes a quantum phase tran- 
sition between different topological phases. The winding 
number Eq. (j^), for this model can be computed nu- 
merically, by discretizing the integral over the BZ. The 
numerical evaluation of as a function of [la and t' is 
shown in Fig. |2|(b), where we set t = i, and A = 2. 
With increasing number of grid points in the BZ the in- 
tegral (H) converges rapidly to an even integer value. In 
this way we obtain the phase diagram shown in Fig. H(a), 
which contains four distinct gapped phases. The non- 
trivial phases {v = ±2) occur if \^s\ < |6i'|. Inclusion of 
finite (but small) second neighbor terms and tg^ shifts 
the phase boundaries, but does not change the topology 
of the phase diagram. 

When and t' are small compared to t [see Fig. 
^(a)], an effective low-energy continuum description can 
be derived by expanding Hamiltonian (^) around the four 
nodal points Ki/2,±- Rescaling momenta as tkz/2 —^ kz, 
t^/3ky/2 ky, and 2^/?t/S.kx — > kx, and performing a 
unitray transformation, we find that the low-energy ex- 
pansion around the node Ki± of the off-diagonal block 
in Eq. (^^ is given by 



D{q) = I3iay [q ■ a ^ i{fis + 6t')j^ 



(9) 



where q = {Qx, ^Qy, ^Qz) denotes the deviation of the 
momentum from Ki±. Eq. (||) is identical to the class CI 
topological Dirac superconductor constructed in Ref . . 
Here, we have introduced the five gamma matrices a^, 
(3, and 7^, which are in the Dirac representation given 



by a. 



Tx, P = Tz, and 7^ = r^r, with Tx^y^z being 



another set of Pauli matrices. Consequently, the winding 
number within the continuum description for the nodes 

Ki,± is 



2 \fis + 6t'\ 



(10) 



where the prefactor 1/2, which is an articfact of the low- 
energy continuum approximation, will be altered once 
information about the structure of Bloch wavefunctions 




FIG. 3: (Color online) (a) Surface Brillouin zone for a slab 
parallel to the (111) surface. Panels (c) and (d) show the 
two-dimensional band structure for the model defined by Eqs. 
(§' terminated by a (111) surface, with the momentum 
being varied along the line indicated in panel (a). For (t', /i) = 
(1.2, 0) [panel (c)] there are surface states which cross the bulk 
energy gap, thereby realizing a surface Dirac fermion state as 
shown in (b). In the trivial phase, for {t',jj,) = (0,4.4) [panel 
(d)], there are no surface states crossing the energy gap. 



at high energy (located away from the Dirac points) is 
supplemented. The low-energy Dirac Hamltonian for the 
nodes i^2,± is related by symmetry to the result for the 
nodes Ki± by simultaneously interchanging kx with ky 
and replacing the mass term fj.s + Gt' with Hs — 6t' , 



D{q) = (3iay [q-a- t(/i, - 6t')7^] 



(11) 



where q = {qy, —Qx, —Qz)- Similarly, the winding number 
within the continuum description for the nodes if2,± is 



1^2 



2 l/Xs 



6t') 



6t' 



X 2. 



(12) 



Interestingly, the winding number obtained from the con- 
tinuum description, ly ~ 1^1 + 1^2, reproduces the phase 
diagram of the lattice model correctly. 

A physical consequence of the non-zero winding num- 
ber ly is the appearance of zero-energy surface Andreev 
bound states with Dirac dispersion. To study these sur- 
face states we solve model (|l|) in a slab geometry, and 
compute the energy bands for a slab parallel to the (111) 
surface, both in the trivial and nontrivial phases. As 
shown in Fig. |[ where we set t — 4 and A ~ 2, there are, 
in addition to the bulk states, surface Dirac states which 
cross the band gap. When the bulk topological invariant 
is ~ 2 [Fig. ||(c)] there are two surface Dirac states, 
whereas there is no such state when the bulk topological 
invariant is = [Fig. ^(d)]. Hence, the total number 
of Dirac states/cones Nf = 2 is consistent with the bulk 
characteristics v ^ 2. Nf = 2 is the minimal number 
of Dirac cones required by class CI symmetries. Note, 
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however, that in any two-dimensional lattice model sat- 
isfying the CI symmetries, the possible number of Dirac 
cones is an integer multiple of '2,Nf because of a no-go 
theorem analogous to the fermion doubling theorem of 
Ref. Here, the fermion doubling is avoided since the 
two-dimensional system is realized as a boundary of a 3D 
bulk. 

The surface Dirac fermion modes cannot be gapped 
by any deformation of the Hamiltonian respecting time- 
reversal and particle-hole symmetry. Indeed, this is so 
because any perturbation respecting the class CI sym- 
metries takes the form of an SU(2) gauge field, which 
perturbs the surface Dirac fermions in the following way 



(13) 



where aJJ^i^'^''^ £ M. The gapless nature of this four- 
component Dirac fermion is stable against arbitrary val- 
ues of the six real parameters a^, i.e., the non-Abelian 
gauge potential shifts the location of the Dirac node, but 
does not lead to a gap. With the inclusion of random- 
ness, small relative to the bulk energy gap, the surface 
Dirac fermion mode realizes the random SU(2) gauge po- 
tential model discussed in Ref. ||l^. The random SU(2) 
gauge potential, known not to be able to localize the 
Dirac fermions, renormalizes to an exactly solved strong 
coupling renormalization group (RG) fixed point at long 
distances. Under the influence of disorder, the (tunnel- 
ing) density of states p(e) changes from a linear depen- 
dence to p(e) ^ |e|^, with the scaling exponent /i = 1/7 
psf at this fixed point. 

Within the non-linear sigma model approach based on 
n fermionic replicas (n —>■ 0), the Anderson localiza- 
tion physics at the surface of the 3D class CI topological 
superconductor in the presence of disorder is described 
by the principal chiral model on Sp(27T.), supplemented 
by a Wess-Zumino-Wittcn (WZW) term, at 'level-one'. 
(When supersymmetric disorder averaging is used, this 
is 11 the WZW model on Osp(2|2) at level [||.) 
We emphasise that the WZW term cannot appear for any 
two-dimensional disordered systems on a lattice, such as 
disordered d^2_y2-wa.ve superconductors on the square 
lattice. Indeed, WZW terms always cancel when con- 
tributions from different cones are added up. Thus, the 
delocalized feature (non-renormalization of the conduc- 
tivity) on a 2D lattice is observed, at best, as a crossover 
behavior when the inter-cone scatterings are weak. 

As it stands, the lattice model we considered is pre- 
sumably not directly connected to any specific system 
occurring in nature; but it may give insight into the prop- 
erties of real materials that exhibit three-dimensional 
topological phases which preserve both time-reversal and 
spin-rotation symmetry. Promising candidates for the CI 
topological quantum state might be found among some 
of the unconventional superconductors of heavy fermion 
systems. There are various ways how the topologically 



protected surface states in such compounds could be de- 
tected experimentally. First of all, because of spin ro- 
tation symmetry, the surface conductivity for the spin 
current is a well-defined quantity, which is unchanged by 
symmetry preserving perturbations, including disorder. 
Secondly, the density of states of the surface states can 
be probed via tunneling experiments. Finally, we argue 
that by breaking TRS locally at the surface, while keep- 
ing the spin-rotation symmetry intact, one can realize 
the so-called 'spm quantum Hall effect (SQHE)'p^ in 
symmetry class C (unrelated to the 'quantum spin Hall 
insulators' mentioned earlier) at the surface of the CI 
topological superconductor. This can be seen as follows. 
Broken TRS (without breaking SU(2) symmetry) allows 
for the appearance of four additional perturbing poten- 
tials in Eq. (13), describing the surface modes. One can 
check that only one of theses additional potentials gives 
rise to a gap, and hence to a (2D) non-zero Chern inte- 
ger n = ±1; this is one half of the allowed value of the 
Chern number in any 2D system exhibiting the SQHE 
Ip^ . This situation is completely analogous to the half- 
integer surface QHE of the 3D Z2 topological insulator 

S 0- 

In the 3D Z2 topological insulator, this topological 
magneto-electric effect can be described by the effective 
field theory whose action is given by (3-1-1) QED supple- 
mented with ^ TT term [|l^, |l^. Since spin is a good 
quantum number in a singlet superconductor, it is possi- 
ble to describe its spin transport in terms of an external 
SU(2) gauge field. Indeed, the effective field theory of 
the SQHE in TRS-breaking 2D singlet superconductors 
is known to be the SU(2) Chern-Simons theory at an in- 
teger level . Then, following the same reasoning as in 
the case of the 3D Z2 topological insulators, we propose 
that the effective field theory describing class CI topo- 
logical superconductor is the (3-|-l)D SU(2) Yang-Mills 
theory argumented with the theta term 



C = 







327r2 



tr (F^uF^x) 



(14) 



where F^,y is the SU(2) field strength and — tt. 
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